clear;
x = {'t-1';'t';'t+1'};
lwid = 2;

x2 = [-1:0.01:4];
x3 = [1;2;3];
y = zeros(length(x2),1);

F1 = [0.000486
0.000362
0.000204];
F2 = [0.000606
0.000390
0.000274];
F3 = [0.000342
0.000028
-0.000021];


tF1 = [3.66
1.99
1.10];
tF2 = [4.47
2.15
1.49];
tF3 = [1.70
0.17
-0.09];

EF1 = 1.64*(F1./tF1);
EF2 = 1.64*(F2./tF2);
EF3 = 1.64*(F3./tF3);


figure(1)
hold on
fill([x3;flipud(x3)],[F1-EF1;flipud(F1+EF1)],[.8 .8 .8],'linestyle','none','FaceAlpha',0.3);
p1 = plot(x3,F1,'k','LineWidth',lwid);
plot(x2,y,'k')
set(gca,'FontSize',13)
ylabel('Coefficient','FontSize',19);
xlim([0.8 3.2])
ylim([-0.0006 0.0009])
xticks([1:1:6])
yticks([-0.0006:0.0002:0.0009])
grid on
xticklabels(x);
ax = gca;
ax.YAxis.Exponent = 0;
ytickformat('%.4f');
hold off
print(gcf,'Figure3a.tiff','-dtiff','-r300'); 

figure(2)
hold on
fill([x3;flipud(x3)],[F2-EF2;flipud(F2+EF2)],[.8 .8 .8],'linestyle','none','FaceAlpha',0.3);
p1 = plot(x3,F2,'k','LineWidth',lwid);
plot(x2,y,'k')
set(gca,'FontSize',13)
ylabel('Coefficient','FontSize',19);
xlim([0.8 3.2])
ylim([-0.0006 0.0009])
xticks([1:1:6])
yticks([-0.0006:0.0002:0.0009])
grid on
xticklabels(x);
ax = gca;
ax.YAxis.Exponent = 0;
ytickformat('%.4f');
hold off
print(gcf,'Figure3b.tiff','-dtiff','-r300'); 

figure(3)
hold on
fill([x3;flipud(x3)],[F3-EF3;flipud(F3+EF3)],[.8 .8 .8],'linestyle','none','FaceAlpha',0.3);
p1 = plot(x3,F3,'k','LineWidth',lwid);
plot(x2,y,'k')
set(gca,'FontSize',13)
ylabel('Coefficient','FontSize',19);
xlim([0.8 3.2])
ylim([-0.0006 0.0009])
xticks([1:1:6])
yticks([-0.0006:0.0002:0.0009])
grid on
xticklabels(x);
ax = gca;
ax.YAxis.Exponent = 0;
ytickformat('%.4f');
hold off
print(gcf,'Figure3c.tiff','-dtiff','-r300'); 